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ABSTRACT 

The measured values of glitch healing parameter, Q, of the Vela pulsar are found to 
be inconsistent with the starquake mechanism of glitch generation in various neutron 
star (NS) models, based upon the parameterized equations of state (EOSs) of dense 
nuclear matter. Since such models correspond to an unrealistic mass range < 0.5Mq 
for the pulsar, if the observational constraints of the fractional moment of inertia of 
the core component (/corc/-^totai < 0.2) which equals the glitch healing parameter, Q, 
in the starquake model, are imposed on these models. However, we show that these 
observational constraints yield a realistic mass range for NS models, corresponding 
to a core given by the stiffest equation of state (EOS), dP/dE = 1 (in geometrized 
units) and the envelope is characterized by the well known EOS of adiabatic polytrope 
dlnP/dlnp = Fi), if the continuity of the adiabatic speed o f sound (v = y^dP/dE)) 
together with pressure (P), energy-density (E), and the two metric parameters {v, 
and A) is assured at the core-envelope boundary of the models and this boundary is 
worked out on the basis of the 'compatibility criterion' for hydrostatic equilibrium. The 
models yield the stable sequence of NS masses in the range, 1.758M0 < M < 2.2Mq, 
corresponding to the values of glitch healing parameter range, < Q < 0.197, for 
a choice of the 'transition density', E^ = 1.342 x lO^^gcm"^, at the core-envelope 
boundary. The maximum stable value of 2.2Mq in this sequence, in fact, corresponds 
to the lowest possible upper bound on NS masses calculated in the literature, on the 
basis of modern EOSs for NS matter. The models yield the surface redshift ~ 0.6913 
and mass M ~ 2.153M0 for the "central" weighted mean value, Q = 0.12 ± 0.07, 
of the glitch healing parameter of the Vela pulsar. This value of mass can increase 
slightly upto M ~ 2.196Mq, whereas the surface redshift can increase upto the value 
zr ~ 0.7568 (which represents an ultra-compact object (UCO; Zfj > 0.73)), if the 
observational constraint of the upper weighted mean value of Q — 0.19 is imposed 
on these models. However, for the lower weighted mean value of Q ~ 0.05, the mass 
and surface redshift can decrease upto the values of M ~ 2.O52M0 and zr ~ 0.6066 
respectively. These results set the lower bound on the energy of a gravitationally 
redshifted radiation in the rather narrow range of 0.291 - 0.302 MeV. The observation 
of the lower bound on the energy of a 7-ray pulse at about 0.30 MeV from the Vela 
pulsar in 1984 is in excellent agreement with this result, provided this energy could be 
interpreted as the energy of a gravitationally redshifted electron-positron annihilation 
radiation from the star's surface. 

Key words: dense matter - equation of state - stars: neutron - stars: pulsars: indi- 
vidual: Vela 



1 INTRODUCTION 

The glitch (sudden increase in the rotational velocity of a 
pulsar) data yield the important information regarding the 
internal structure of neutron stars (NSs) , since it provide the 
best tool for testing various glitch models which are actually 
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based upon the internal structure of NSs. At present, there 
are two well accepted glitch models available in the litera- 
ture: (i) the starquake (Baym & Pines 1971; Ruderman 1976; 
Alpar et al 1996) and (ii) the vortex unpinning (Anderson & 
Itoh 1975; Alpar et al 1993) models. And, both of them, in 
fact, agreed, in general, upon the same conventional notion 
that a NS may be considered as a two-component structure, 
a superfluid interior core which contains most of the NS's 
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mass surrounded by a rigid crust which contains only a few 
percent of the total mass (here the term 'crust' is used for 
the solid crust plus other interior part of the star (strongly 
coupled to it) right upto the superfluid core; we shall call this 
portion as 'envelope' of the star). The envelope of conven- 
tional NS models is characterized by different equations of 
state (EOSs) in am appropriate sequence below a 'fiduciary' 
transition density, Ei,, and the core is sometimes character- 
ized by the extreme causal EOS, dP/dE — 1 beyond Et, 
in order to calculate an upper bound on NS masses. Each 
member of the stable mass-radius sequence of such models 
may represent a 'realistic' case. 

At present, the most extensive and accurate glitch data 
for Crab and Vela pulsars are available in the literature. 
Crawford & Dcmianski (2003) have collected the all mea- 
sured values of glitch healing parameter, Q, for Crab and 
Vela pulsars, which is defined in starquake glitch model as 
the fractional moment of inertia, i.e. the ratio of the mo- 
ment of inertia of the superfluid core, /core, to the moment 
of inertia of the entire configuration, /total, as 



Q = 



/total 



(1) 



They have calculated that 21 measured values of Q for Crab 
glitches yields a weighted mean of Q = 0.72 ± 0.05, and the 
range of Q > 0.7 encompasses the observed distribution for 
the Crab pulsar. In order to test the starquake model for 
Crab pulsar, they have computed Q (as given by Eq.(l)) 
values for seven representative EOSs of dense nuclear mat- 
ter, covering a range of neutron star masses. Their study 
shows that the much larger values of Q(> 0.7) for the Crab 
pulsar is fulfilled by all, but the six EOSs considered in the 
study corresponding to a 'realistic' neutron star mass range 
1.4 ± 0.2Mo. On the other hand, a weighted mean value of 
the 11 measurements for Vela yields a much smaller value of 
Q{= 0.12±0.07) and the all estimates for Vela agree with the 
likely range of Q < 0.2. Thus, their results are found to be 
consistent with the starquake model predictions for the Crab 
pulsar. They have also concluded that the much smaller val- 
ues of Q < 0.2 for Vela pulsar are inconsistent with the star- 
quake model predictions, since the implied Vela mass based 
upon their models corresponds to a value < O.O5M0, which 
is too low as compared to the 'realistic' NS mass range. 

However, it seems really surprising that if the NS's in- 
ternal structure is described by the same conventional NS 
models (as mentioned above), why different kinds of glitch 
mechanisms are required for the explanation of a glitch! It 
is apparent from the study of weighted mean value of the 
glitch healing parameter, Q, for the Vela pulsar that one 
would require a two-component model of NS such that more 
than eighty percent of the total moment of inertia of the 
star should remain confined in the envelope region in order 
to have the starquake model explanation of Vela glitches. 
The present study deals with the construction of such a 
model which is possible if wo impose the conditions: (i) to- 
gether with other variables (pressure, energy-density, both 
of the metric parameters v and A), the (adiabatic) speed 
of sound v{= ^/IdP/dE)) should also remain continuous at 
the corc-cnvclopo boundary of the model governed by an 
EOS of 'adiabatic' polytrope, dlnPjdlnp = Fi, in the enve- 
lope region and the extreme causal EOS, {dP/dE) = 1, in 
the core region respectively, (ii) The models obtained in (i) 



should follow the 'compatibility criterion' that for every as- 
signed value of the ratio of central pressure to central energy- 
density (cr = Po/Eo), the compactness ratio w(= M/R; to- 
tal mass to radius ratio of the static configuration in ge- 
omctrizod units) of the models should remain less than or 
equal to the compactness ratio of the corresponding sphere 
of homogeneous density distribution, in order to assure the 
condition of hydrostatic equilibrium (Negi & Durgapal 2001; 
Negi 2004a). 

The reason for assigning the above mentioned gamma- 
law EOS for densities below the fiduciary transition density 
Eh, is not only because of the fact that matching of speed of 
sound (together with other variables) at the core-envelope 
boundary is possible in an analytical maimer (Negi & Dur- 
gapal 2000), but also because such models (for various as- 
signed values of constant Fi) provide the upper bound on 
NS masses independent of the EOS of the envelope, if the 
condition of 'compatibility criterion' is included (Negi 2005). 
In the present context, however, the value of Fi follows itself 
from the matching conditions at the core-envelope boundary 
on the basis of 'compatibility criterion' and may be looked 
upon as an 'average' (constant) value of Fi below the density 
range, Eb, if this range could have been specified by various 
EOSs like, WFF (Wiringa, Fiks & Fabrocini 1988), FPS 
((Lorenz, Ravenhall & Pethick, 1993), NV (Negele & Vau- 
therin 1973), or BPS (Baym, Pethick & Sutherland 1971) in 
an appropriate sequence, as are frequently used by various 
authors in the conventional models of NSs (sec, e.g. Kalogera 
& Baym 1996; Friedman & Ipser 1987), since the choice of 
the fiduciary transition density, E^ = 1.342 x lO^^gcm"^ 
(which is higher than 4Enm, where Enm = 2.7 x lO^'' gcm~^ 
denotes the nuclear matter saturation density) yields an up- 
per bound on stable NS masses, Mmax = 2.2Mq, for our 
model which, in fact, corresponds to the lowest possible up- 
per bound on NS masses, independent of Eb if Eb > AEnm, 
calculated by Kalogera & Baym (1996) on the basis of mod- 
ern EOSs for NS matter, fitted to experimental nucleon- 
nucleon scattering data and the properties of light nuclei. 
This result also shows consistency with the reasoning men- 
tioned in the previous sentence regarding the imposition of 
the 'compatibility criterion' on NS models. The reproduction 
of the measured values of the glitch healing parameter, Q, for 
the Vela pulsar on the basis of the present study indicates 
that this pulsar should be much compact (approximately 
twice) as compared to that of the Crab (see, e.g. Negi 2005) 
and suggests a further implication (discussed under section 
4), if starquake is considered to be a viable mechanism for 
glitch generation in all pulsars. 



2 METHODOLOGY 

The metric for spherically symmetric and static configura- 
tions can be written in the following form 



ds 



e dt 



A , 2 

e ar 



(2) 



where v and A are functions of r alone. The units are as- 
signed in such a manner that the Newtonian gravitation 
constant (G) and the speed of light in vacuum (c) both 
become unity. The Oppenheimer-Volkoff (O-V) equations 
(Oppenheimer & VolkofF 1939), resulting from the Einstein's 
field equations for the systems with isotropic pressure P and 
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energy-density E can be written as 



P' 



-{P + E)[4TTPr'^ + m]/r(r - 2m), 



^'/2 = -P'/{P + E), 
m'{r) = AnEr'^ ; 



(3) 
(4) 
(5) 



where the prime denotes radial derivative and m(r) is the 
mass contained within the radius r 



m(r) 



-f 

Jo 



4nEr dr. 



The core (0 < r < b) of the present model can be writ- 
ten in the following form 



P = {E-Es 



(6) 



where Es is the surface density of the configuration char- 
acterized by the extreme causal EOS. While the envelope 
(6 < r < R) is given by the EOS 



P = Kp'' 



(7) 



{E-p) = P/{r,-l). 



where K is a constant to be worked out by the matching 
of various variables at the core-envelope boundary and p 
and Fi represent respectively, the rest-mass density and the 
(constant) adiabatic index (sec, e.g., Tooper 1965). 

At the core-envelope boundary, r = b, the continuity of 
P{= Pb),E{= Eb), and r{= rb) require ( Negi & Durgapal 
2000) 



K = Pb/{Eb - [Pb/{Ti 
where Fi is given by 



(8) 



Fi = [{P + E)/P]{dP/dE). 

The continuity of (dP/dE), at the core-envelope bound- 
ary requires 



Fi = 1 + (Eb/Pb). 



(9) 



Thus, the contirmity of (dP/dE), together with other vari- 
ables (P, E, v, and A) is ensured at the core-envelope bound- 
ary of the static and spherically symmetric configuration. 

The coupled Eqs.(3), (4), (5), are solved for the model 
by considering Eq.(6) in the core and Eq.(7) in the envelope 
alongwith the boundary conditions (8) and (9) at the core- 
envelope boundary, r = b, and the boundary conditions, P — 
E = ,m{r = R) = M, e" = e^^ = (1 - 2M/i?) = (1 - 2u) 
aX r = R, at the surface of the configuration, such that for 
each possible value of cr, the compactness parameter of the 
whole configuration always turns out to be less than or equal 
to the compactness parameter of the corresponding sphere 
(with the same cr) of the homogeneous density distribution. 
It is seen that this condition is fulfilled if the minimum value 
of the ratio of pressure to energy-density, Ph/Eb, at the core- 
envelope boundary reaches about 2.92 x 10~^. The results 
of the calculations are presented in Table 1 and 2 respec- 
tively, while the mass-radius diagram is shown in Fig.l for 
a choice^ of Eb — 1.342 x 10^^ gcm~^ (which is higher than 
AEnm, where E„m = 2.7 x lO'^* gcm~^ is the nuclear matter 

^ this choice of transition density and its subsequent implications 
are discussed under section 4 of the present study. 



saturation density) . We find that the first maxima in mass is 
reached to a value of 2.2M0 for a o value about 0.6285 which 
is the evidence that the models become pulsationally stable 
up to the maximum value of mass Mmax = 2.2M0. The 
radius corresponding to this maximum mass is obtained as 
9.587 km. Thus, the maximum compactness, Umax, for this 
stable model yields the value ~ 0.3389, as shown in Table 
1. This upper bound is found to be fully consistent with the 
exact absolute upper bound on compactness ratio of NSs 
compatible with causality and pulsationally stability (sec, 
e.g. Negi 2004b), and provides evidence regarding the ap- 
propriateness of the model with the Fi = constant envelope. 
The binding energy per unit rest-mass a,- [= {Mr — M)/Mr\ 
where Mr is the rest-mass (sec, e.g. Zcldovicli & Novikov 
1978)] also approaches its first maxima of about 0.2314 for 
the first maximum value of mass up to which the configura- 
tions remain pulsationally stable. 



3 AN APPLICATION OF THE MODELS TO 
VELA PULSAR 

For slowly rotating configurations like Vela pulsar (rotation 
velocity, fi, about 70 rad sec"'^) the moment of inertia may 
be calculated in the first order approximation that appears 
in the form of Lense-Thirring frame dragging-effect. In such 
situations, however, it appears very useful to use an approx- 
imate, but very precise empirical formula which is based on 
the numerical results obtained for thirty theoretical EOSs 
of dense nuclear matter. For NSs, the formula yields in the 
following form (Bejger & Haensel 2002) 

, 2 



-{l + hx)MR\ x>Q.l 



(10) 



where x is the compactness parameter measured in units 
of [M0(km)/km], i.e. 



M/R 



1.477' 



(11) 



Me/km 

Eq.(lO) is used, together with coupled Eqs.(3 - 5), to 
calculate the fractional moment of inertia given by Eq.(l) 
and the moment of inertia of the entire configuration for 
models presented in Table 1 and Fig.l respectively. The re- 
sults are shown in Table 2 and Fig. 2 respectively. For the 
central weighted mean value Q ~ 0.12 for the Vela pulsar, 
Fig.2 yields the mass value M ~ 2.153M0, radius R ~ 9.780 
km, and surface redshift Zr ~ 0.6913 respectively. The cor- 
responding core mass, Mb, turns out to be about O.SIOMq 
as shown in Table 2 which is about 37.6% of the total mass of 
the entire configuration. This value of Vela mass can exceed 
upto a value about 2.196M0 if the upper limit of the central 
weighted mean value of Q ~ 0.19 is considered. The corre- 
sponding values of surface redshift and core mass can exceed 
upto the values 0.7568 and 1.O46M0 respectively. This value 
of surface redshift {zb. — 0.7568), in fact, represents an ultra- 
compact object (UCO; zjj > 0.73) which are entities of inter- 
est (see, e.g. Negi & Durgapal 1999; and references therein). 
For the lower limit of the 'central' weighted mean value of 
Q ~ 0.05, the Vela mass can reduce upto a value of about 
2.O52M0. The corresponding values of surface redshift and 
core mass turn out to be 0.6066 and O.489M0 respectively. 
Thus, the core mass of the structure varies between 23.8 % 
to 47.6% of the total mass, for the lower limit of Q ~ 0.05 
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Figure 1. Mass- Radius diagram of the models as discussed in 
the text for an assigned value of the transition density E = Ef, = 
1.342x10"'''' g cm^'^ at the core-cnvelopc boundary. The minimum 
value of the ratio of pressure to energy-density, (Pi,/Ei,), at the 
core envelope boundary is obtained as 2.92 X 10^^, sueh that for 
an assigned value of a, the inequality u < ujj is always satisfied, 
as shown in Table 1. 




Figure 2. Mass (M/Mq) vs. fractional moment of inertia Q{= 
^core/Aotai) for the models as discussed in the text and presented 
by Table 1 and 2 respectively. 



to the upper limit of Q ~ 0.19 respectively. For Q values 
larger than about 0.197 (i.e., the point of maximum mass 
Mmax = 2.2Mq), the structures become pulsationally un- 
stable, whereas the minimum stable Vela mass corresponds 
to a value about 1.758 Mq as Q — > 0. 



4 DISCUSSION 

In the previous study, we considered the values of con- 
stant Fi = (4/3), (5/3), and 2 respectively for the density 
range below the fiduciary transition density Ei, on the basis 
of 'compatibility criterion' in order to construct the star- 
quake models for Crab pulsar (Negi 2005). If (i) the ob- 
servational constraint of the glitch healing parameter, and 
(ii) the observational constraint of the recently evaluated 
value of the moment of inertia for the Crab pulsar were 
combined together with the 'compatibility criterion' men- 
tioned above, the model with Fi = (5/3) envelope itself 
yielded the value of transitions density, ~ 2.7 x lO'^* gcm~^, 
the nuclear matter saturation density at the core-envelope 
boundary. This value of Ei, yields the upper bound on NS 
masses Mmax — 4.IM0, independent of the EOS of the en- 
velope (Negi 2005). Not only this value is found fully consis- 
tent with the upper bound on NS masses obtained by using 
modern EOSs for NS matter, fitted to experimental nucleon- 
rmcloon scattering data and the properties of light nuclei, 
consistent with causality and dynamical stability (see. Fig. 2 
of Kalogera & Baym 1996, for a fiduciary transition density 
Eb = 2.7 X lO^^gcm"*), the upper bound on compactness 
ratio of these models {umax — 0.34) was also found fully 
consistent with the 'absolute' upper bound on compactness 
ratio of NSs consistent with causality and pulsational stabil- 
ity (Negi 2004b) . The results of this study are also supported 
by observations, since the existed value of the upper bound 
on the energy of a gravitationally redshiftod 7-ray line at 
about 0.40 MeV from the Crab pulsar is found consistent 
with the predicted value in the energy range of about 0.414 
- 418 MeV (see, Negi 2005; and references therein). 

In the previous study, the value of the ratio of pres- 
sure to density at the core-envelope boundary, Pt/Et, on 
the basis of 'compatibility criterion' was calculated as ~ 
1.065 X 10~^, which represents the region very close to the 
surface of the star (thin crust), and does not require the 
necessary continuity of the speed of sound (together with 
other variables) at this boundary which (obviously) corre- 
sponds to a region of relatively low density. Since, we were 
concerned with the larger values of the glitch healing pa- 
rameter, Q > 0.70, for the Crab pulsar which demanded 
that more than 70% of the moment of inertia should re- 
main attached to the core component. However, the much 
smaller values of the glitch healing parameter for Vela pul- 
sar, Q < 0.20, require that more than 80% of the moment of 
inertia should remain confined to the envelope region which 
obviously require a much larger value of Pb/Ef, (thick crust) 
as compared to the previous case, i.e. a much deeper region 
of higher density inside the crust where one can naturally 
assume that the speed of sound approaches the speed of 
light and is contirmous at some 'fiduciary' density at the 
boundary of extreme causal EOS in the core. In the present 
context of the Vela pulsar, we do not have any informar 
tion about the Vela's moment of inertia, as we did have in 
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the case for the Crab. However, if the Vela and the Crab 
both are members of the same conventional sequence of the 
NS models which terminates at the maximum value of NS 
meiss, then we can assign the value of 4.1Mq as an upper 
bound of NS masses obtained in the previous study, in order 
to constrain the present models. This assignment gives the 
Vela's mass about 4.01Mq for the 'central' weighted mean 
value of the glitch healing parameter Q ~ 0.12 which looks 
too large, since none of the observational studies claim such 
higher values of mass for the Vela pulsar. The correspond- 
ing value of transition density, turns out to be about 
3.864 X 10^*gcm~^, which also seems unlikely for the nor- 
malization of the Pb/Ef, ratio of ~ 2.92 x 10~^ calculated for 
the present models. Notwithstanding, the result of the pre- 
vious study which is found fully consistent with the present 
one is the upper bound on NS compactness, u ~ 0.34. 

This finding reveals that if starquake is considered to 
be a viable model for glitch generation in the Vela pulsar^ 
together with the Crab (and, in fact, in all pulsars!) then 
it would lead to the following implication that there might 
exist a parallel sequence of the stable NS masses composed 
of Vela-like pulsars, together with the conventional sequence 
which composes the Crab- like pulsars, such that this parallel 
sequence begins with somewhat higher value of the compact- 
ness ratio as compared to the conventional sequence but ter- 
minates with almost the same value of maximum compact- 
ness (e.g., the parallel sequence considered in the present 
study begins with a value of compactness ratio about 0.2628 
which is higher than the corresponding value of compactness 
ratio of conventional sequence undertaken by Negi (2005), 
but both of these sequences terminate with almost the same 
value of maximum compactness, u ~ 0.34) consistent with 
the absolute upper bound on NS compactness (Negi 2004b). 
However, such a sequence would correspond to different val- 
ues of minimum and maximum stable masses as compared to 
the conventional sequence. Apparently, the minimum mass 
of this sequence would be higher than the minimum mass 
of the conventional sequence, whereas the maximum mass 
of this sequence would correspond to a value less than the 
maximum mass of the conventional sequence. 

Kalogera & Baym (1996) have shown that the lowest 
possible upper bound on NS masses corresponds to a value 
of 2.2Mq independent of Et, provided the fiduciary tran- 
sition density, Ei,, satisfies the condition Et, > AEnm- Our 
models, in fact, satisfy this condition, since for the lowest 
possible upper bound of 2.2Mq, we obtain the value of tran- 
sition density, Et = 1.342 x 10^^ g cm^"^, which is higher than 
4-Enm,. Apparently, this value of transition density yields 
the 'minimum' possible mass for Vela pulsar on the basis of 
present models as 2.153Mo for the 'central' weighted mean 
value of Q ~ 0.12. This is discussed earlier under sections 
2 and 3 respectively. Thus, in view of the maximum mass 
and density regime, the results of the present study are fully 

The other problems associated with the starquake model for 
Vela pulsar, like the large size of glitehes {AQ/Q ~ 10"'^), and 
the observations of significant amount of change in X-ray flux soon 
after the glitch occurs (see, e.g. Crawford & Demiahski (2003); 
and references therein ) are not discussed in the present paper. 
However, in view of the higher values of compactness ratio of the 
present models, the future study in this regard may provide some 
explanation. 



consistent with the model based on the modern theory of nu- 
clear interactions and underline the appropriateness of the 
use of gamma-law EOS in the envelope, the necessity of the 
continuity of speed of sound across the boundary, and the 
imposition of the 'compatibility criterion', in order to have 
a starquake model explanation for Vela-like pulsars. 



5 RESULTS AND CONCLUSIONS 

The study shows that if starquake is responsible for the 
glitch generation in Vela pulsar, the surface redshift, zr, 
for NS model of Vela pulsar should correspond to a value 
of 0.6912, if the observational constraint of the 'central' 
weighted mean value of the glitch healing parameter Q ~ 
0.12 is imposed. This value gives the 'minimum' mass of 
the Vela pulsar M ~ 2.153M0 for the 'transition density' 
Eb = 1.342 X 10^® gcm~^ at the core-envelope boundary, if 
the lowest possible upper bound on NS masses is considered 
to be 2.2Mq . The mass is slightly increased upto the value 
M ~ 2.196M0, but the surface redshift can increase upto 
the value zr ~ 0.7568 if the observational constraint of the 
upper weighted mean value of the glitch healing parameter 
Q ~ 0.19 is imposed. However, for the lower weighted mean 
value of the glitch healing parameter Q ~ 0.05, the mass 
and surface redshift can decrease upto zr ~ 0.6066 and 
M ~ 2.O52M0 respectively. For the upper weighted mean 
value of Q ~ 0.19, the structures represent ultra-compact 
objects (UCO) which are entities of important astrophysi- 
cal interest (see, e.g. Negi & Durgapal 1999; and references 
therein). In the limiting case of Q — > 0, the minimum mass 
for Vela corresponds to the value of 1.758M0. The confir- 
mation of these results require an observation of the lower 
bound on the energy of a gravitationally redshifted 7-ray 
line in the narrow energy range about 0.291 - 0.302 MeV 
from the Vela pulsar which is in excellent agreement with 
the observation of the lower bound on the energy of a 7- 
ray pulse of 0.30 MeV from the Vela pulsar (Tumcr ct al. 
1984), if this energy could be interpreted as the energy of a 
gravitationally redshifted electron-positron annihilation rar 
diation from the star's surface. 
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Table 1. Various values of mass (M/Mq), radius _R(krn), compactness ratio (u), binding-energy per particle (ctr), corresponding to the 
models discussed in the text for different assigned values of a. The compactness ratio for homogeneous density distribution is represented 
by Ufi- The slanted values correspond to the limiting case upto which the configuration remains pulsationally stable. 



(Po/Eo) M/Mq R{km) ar u 



0, 


.29202 


1, 


.757659 


9. 


.877204 


0. 


.169796 


0. 


.262834 


0. 


.262854 


0. 


.29800 


1, 


.772991 


9. 


.884901 


0. 


.171393 


0. 


.264920 


0. 


.265166 


0, 


.35452 


1, 


.903785 


9. 


.910778 


0, 


.186960 


0. 


.283720 


0. 


.284570 


0, 


.39993 


1, 


.991678 


9. 


,920844 


0, 


.198449 


0, 


.296518 


0, 


.297502 


0, 


.43923 


2, 


.052262 


9. 


,896808 


0, 


.206909 


0, 


.306280 


0, 


.307194 


0, 


.45699 


2, 


.077665 


9. 


,892000 


0, 


.210533 


0, 


.310222 


0, 


.311190 


0, 


.49397 


2 


.120681 


9, 


,854696 


0, 


.217014 


0, 


.317843 


0, 


.318832 


0, 


.53331 


2 


.153382 


9, 


,780031 


0, 


.222426 


0, 


.325208 


0, 


.326097 


0, 


.54099 


2, 


.159962 


9. 


.772185 


0. 


.223498 


0. 


.326464 


0. 


.327424 


0, 


.55999 


2, 


.171257 


9. 


.729357 


0. 


.225572 


0. 


.329615 


0. 


.330584 


0. 


.59146 


2, 


.189320 


9. 


.678408 


0. 


.228886 


0. 


.334107 


0. 


.335476 


0, 


.61902 


2, 


.196170 


9. 


.597157 


0, 


.230666 


0. 


.337990 


0. 


.339441 


0. 


.62850 


2 


.200000 


9. 


.587457 


0, 


.231448 


0, 


.338922 


0, 


.340741 


0, 


.62936 


2, 


.197304 


9. 


,563581 


0, 


.231138 


0, 


.339352 


0, 


.340858 


0, 


.63148 


2, 


.198152 


9. 


,560925 


0, 


.231311 


0, 


.339577 


0, 


.341144 


0, 


.65148 


2, 


.200457 


9. 


,506329 


0, 


.232207 


0, 


.341885 


0, 


.343769 


0, 


.69998 


2, 


.195136 


9. 


,351309 


0, 


.232781 


0, 


.346713 


0, 


.349633 


0, 


.75294 


2, 


.171828 


9. 


.151000 


0. 


.230592 


0. 


.350540 


0. 


.355328 



Table 2. Various values of mass {M/Mq), radius iZ(km), core mass {Mi,/Mq), core radius i?i,(km), and fractional moment of inertia 
^core /Aotal ) corresponding to the models discussed in the text for different assigned values of a. The corresponding surface and central 
redshifts are represented by zn and zq respectively. The slanted values correspond to the limiting case upto which the configuration 
remains pulsationally stable. 



(./i)/-t") 


M/.\.l . 


./i'(kni) 


..W,,/-\./ . 


./i'/,(km) 


./,-, 


.r,.//,,.,.l 








-0 


0. 


,29202 


1, 


.757659 


9, 


.877204 


0, 


.000001 


0, 


.084931 


0, 


.000000 


0, 


.451973 


0, 


.961475 


0. 


,29800 


1, 


.772991 


9, 


.884901 


0, 


.006151 


1, 


.301636 


0, 


.000032 


0, 


.458402 


0, 


.983101 


0. 


,35452 


1, 


.903785 


9, 


.910778 


0, 


.171574 


3, 


.889230 


0, 


.008641 


0, 


.520468 


1, 


.202564 


0. 


.39993 


1. 


.991678 


9. 


.920844 


0, 


.341592 


4, 


.844349 


0. 


.027604 


0, 


.567552 


1, 


.394268 


0. 


.43923 


2. 


.052262 


9. 


.896808 


0, 


.489306 


5. 


.412447 


0. 


.050835 


0, 


.606561 


1, 


.573780 


0. 


.45699 


2. 


.077665 


9. 


.892000 





.554063 


5, 


.617970 


0. 


.062644 


0, 


.623161 


1, 


.658789 


0. 


.49397 


2, 


.120681 


9. 


.854696 





.683071 


5, 


.969921 


0. 


.089515 


0, 


.656770 


1, 


.846732 


0. 


.53331 


2, 


.153382 


9, 


.780031 





.809999 


6, 


.255151 


0. 


.120661 


0, 


.691314 


2, 


.065669 


0. 


.54099 


2, 


.159962 


9, 


.772185 


0, 


.833445 


6, 


.301950 


0, 


.126634 


0, 


.697422 


2, 


.110140 


0. 


,55999 


2, 


.171257 


9, 


.729357 


0, 


.889491 


6, 


.406664 


0, 


.142236 


0, 


.713049 


2, 


.225541 


0. 


.59146 


2, 


.189320 


9, 


.678408 


0, 


.976141 


6, 


.548950 


0, 


.167190 


0, 


.736086 


2, 


.426520 


0. 


.61902 


2, 


.196170 


9, 


.597157 


1, 


.045560 


6, 


.644820 


0, 


.190181 


0, 


.756766 


2, 


.621490 


0. 


.62850 


2, 


.200000 


9, 


.587457 


1, 


.068060 


6 


.672101 


0. 


.197133 


0, 


.761842 


2, 


.688782 


0. 


.62936 


2. 


.197304 


9. 


.563581 


1, 


.070062 


6. 


.674508 


0. 


.198876 


0, 


.764196 


2, 


.698936 


0. 


.63148 


2. 


.198152 


9. 


.560925 


1, 


.074988 


6, 


.680161 


0. 


.200427 


0. 


.765435 


2, 


.714632 


0. 


.65148 


2. 


.200457 


9. 


.506329 


1, 


.119576 


6. 


.727057 


0. 


.216372 


0. 


.778276 


2, 


.871016 


0. 


.69998 


2. 


.195136 


9. 


.351309 


1, 


.214118 


6. 


.793253 


0. 


.254274 


0, 


.806058 


3, 


.299069 


0. 


.75294 


2. 


.171828 


9. 


.151000 


1, 


.294603 


6. 


.792263 


0. 


.293308 


0. 


.829036 


3, 


.871991 



